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ABSTRACT: Knowledge of the frontier levels’ alignment 
prior to photo-irradiation is necessary to achieve a com¬ 
plete quantitative description of H 2 O photocatalysis on 
TiO 2 ( 110 ). Although H 2 O on rutile TiO 2 ( 110 ) has been 
thoroughly studied both experimentally and theoretically, a 
quantitative value for the energy of the highest H 2 O occu¬ 
pied levels is still lacking. For experiment, this is due to the 
H 2 O levels being obscured by hybridization with TiO 2 ( 110 ) 
levels in the difference spectra obtained via ultraviolet pho¬ 
toemission spectroscopy (UPS). For theory, this is due to in¬ 
herent difficulties in properly describing many-body effects 
at the H 20 -Ti 02 (l 10) interface. Using the projected density 
of states (DOS) from state-of-the-art quasiparticle (QP) GqWq, we disentangle the adsorbate and surface contributions to the complex 
UPS spectra of H 2 O on TiO 2 ( 110 ). We perform this separation as a function of H 2 O coverage and dissociation on stoichiometric 
and reduced surfaces. Due to hybridization with the TiO 2 ( 110 ) surface, the H 2 O 3ai and Ibi levels are broadened into several peaks 
between 5 and 1 eV below the Ti 02 (l 10) valence band maximum (VBM). These peaks have both intermolecular and interfacial bonding 
and antibonding character. We find the highest occupied levels of H 2 O adsorbed intact and dissociated on stoichiometric TiO 2 ( 110 ) 
are 1.1 and 0.9 eV below the VBM. We also find a similar energy of 1.1 eV for the highest occupied levels of H 2 O when adsorbed 
dissociatively on a bridging O vacancy of the reduced surface. In both cases, these energies are significantly higher (by 0.6 to 2.6 eV) 
than those estimated from UPS difference spectra, which are inconclusive in this energy region. Finally, we apply self-consistent 
QPGW (scQPGWl) to obtain the ionization potential of the H 2 O-TiO 2 ( 110 ) interface. 



1. INTRODUCTION 

The photooxidation activity of a surface is determined by the 
interfacial level alignment between the occupied adsorbate lev¬ 
els and those of the substrate.^ Water photooxidation on TiO? 
has attracted enormous attention^^®^ for energy applications!^^® 
based on H 2 production.® This reaction also plays an important 
role in photocatalv tic env ironmental remediation and surface self¬ 
cleaning/sterilizing.®® This is because the resulting hydroxyl 
radicals are the key intermediates in the oxidative degradation of 
organic species.®^ To understand water photooxidation, it is nec¬ 
essary to understand the interfacial level alignment between the oc¬ 
cupied levels of H 2 O and the Ti 02 substrate.® 

Experimentally, the most common approach to access the adsor¬ 
bate levels is to take the difference between the covered and clean 
surface spectra from photoemission spectroscopy. However, when 
the adsorbate and surface levels are strongly hybridized, it becomes 
difficult to disentangle the adsorbate and surface contributions to 
the UPS spectra using only the difference spectra.® For example. 


shifting of the surface levels due to hybridization or band bending 
may completely obscure the adsorbate levels.' — 'Further, the adsor¬ 
bate levels near the valence band maximum (VBM) are the most 
likely to be obscured. It is precisely these levels that are most im¬ 
portant for photooxidation processes. Using a theoretical approach, 
one can directly disentangle the molecular levels by projecting the 
density of states (DOS) of the interface onto the atomic orbitals of 
the molecule. Altogether, this makes a robust theoretical approach 
necessary to accurately predict the alignment of the adsorbate and 
substrate levels, and separate the adsorbate and surface spectra. 

A robust theoretical treatment requires quasiparticle (QP) GqWq 
to capture the anisotropic screening of the electron-electron in¬ 
teraction at the interface.®® As previously demonstrated for 
CH 3 OH on Ti 02 (l 10), QP GqWq is necessary to obtain even a qual- 
itative description of the level alignment.' " 'For this interface, the 
occupied levels of the molecule are only weakly hybridized with 
the surface levels. This allowed an unambiguous comparison to the 
photoemission difference spectrum.' — 'However, for H 2 O on rutile 
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TiO2(110), this is not the case. 

The occupied molecular levels of H2O on single crystal ru¬ 
tile TiO2(110) hav e been probed via ultraviolet photoemission 
spectroscopy (UPS)®^^^ and metastable impact electron spec- 
trosconv tMlRSV^ These exneriments were nerformed under iil- 


troscopy (MIES). ^ These experiments were performed under ul- 
trahi^ vacuum (UHV) conditions from low to room tempera¬ 
ture “from 0.01 to 100 L H2O exposure,^ and for various sur¬ 
face preparations resulting in either reduced TiO2-jcQ10) with sur- 
face oxygen defects or “nearly-perfect” TiO2(110).' — 'Altogether, 
these experiments have addressed the long-standing controversy as 
to where and how H2O adsorbs and dissociates on Ti02(l 10).!^^^^ 

At 150 K the photoemission difference spectrum between H2O 
covered and clean TiO2(110) surfaces consists of three peaks, 

which are attributed to intact H2O adsorbed on Ti coordinately 

i2n 

unsaturated sites (Ticus)- —' Upon heating to 300 K, the difference 
spectrum’s three-peak structure evolves into a two-peak structure, 
which is attributed to dissociated H2O adsorbed on bridging O va¬ 
cancies i.e., Ohrii surface species.^ This assignment of the 

UPS spectra to intact (I) ll 20 @Ticus or dissociated (D) H20@0^^^ 
is based on the peak energy separations being consistent with those 
reported for H2cJ^in gas phase or OH“ in NaOH.^ 

A comparison to the H2O and OH“ peaks is robust for the molec¬ 
ular levels that lie below and have little hybridization with the sur¬ 
face DOS. However, the adsorbate levels that lie within the surface 
valence band may significantly hybridize with the surface, with a 
single molecular level contributing to many interfacial levels. These 
interfacial levels are thus not easily associated with H2O and OH“ 
levels. This is exacerbated by the mixing of the molecular levels 
due to symmetry breaking at the interface. As a result, “between 5 
and 8 eV” below the Fermi level, experimentally they “are unable 
to produce reliable difference structures” from the UPS spectra ob¬ 
tained for “nearly-perfect” TiO2(110) exposed to H2O at 160 K.^ 

Using the QP GqWq H2O projected DOS (PDOS), we have dis¬ 
entangled the adsorbate and surface contributions to the UPS spec¬ 
tra within this difficult energy range. This has been done as a func¬ 
tion of H2O coverage and dissociation on stoichiometric and re¬ 
duced surfaces. In so doing, we provide quantitative values for 
the energies of the highest H2O occupied levels, prio r to p hoto- 
irradiation, for a number of experimentally relevan t H2O- 
TiO2(110) structures. 

To directly compare to red-ox potentials, the important quantities 
for determining photoelectrocatalytic a ctivity , one needs the align¬ 


ment relative to the vacuum level, With this, one obtains 

the ionization potential directly from + Evac- To obtain a 

more accurate absolute level alignment, we employ our recently in¬ 


troduced self-consistent QP G 


lignmen 

vpEl 


technique scQPGWl.' 
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The presentation of the results is organized as follows. First, we 
focus on the H2O levels that lie below and have little hybridiza¬ 
tion wit h the substrate DOS. This is done for intact^ 0 @TicM^ in 
Section ' 


Section 


3.1 


3.3 


3.2 


Further, in 


and dissociated H20@0^^ in Section 
we shown that these results are rather independent of 
the choice of xc-functional. In so doing we provide evidence for a 
robust semi-quantitative agreement with the UPS difference spectra 
for the adsorbate levels for which an unambiguous comparison with 
the experiment is possible. For a more complete understanding of 
the UPS experiments, in Section [T 4 | we analyze the H2O PDOS for 
a variety of other H2O structures on the stoichiometric and reduced 
surfaces. These may form under different experimental conditions 
and surface preparations. In Section [T 5 ] we focus on the highest 
H2O occupied levels, which are significantly hybridized with the 
substrate DOS. The success of the QP GqWq PDOS strategy for the 
lower-energy part of the UPS difference spectra provides support 
for our results in this difficult spectral region, where a straightfor¬ 
ward comparison with experiment is not possible. Finally, in Sec¬ 


tion | 3 . 6 | we employ scQPGWl to obtain an improved absolute level 
alignment relative to Evac, and thus estimate the ionization potential 
of the H2O-TiO2(110) interface. 

2. METHODOLOGY 

Our QP Go Wo calculations®®^ have been performed using vasp 
within the projector augmented wave (PAW) scheme.®^ The Go Wo 
calculations are based on Kohn-Sham wavefunctions and eigenen- 
ergies from density functional theory (DFT) obtained using a gen¬ 
eralized gradient approximation (PBE)®^ for the exchange corre¬ 
lation (xc) -functional.® The dependence of the QP Go Wo DOS 
and PDOS on the DFT xc-functional has been tested for 1 ML 
intact H20@TicM5 of stoichiometric TiO2(110) and ViML dissoci¬ 
ated H20@0^f of defective Ti02-i/4(l 10 ) with ViML of O^f. For 
these structures, GoWn calculations based on the local density ap¬ 
proximation (LDA),® van der Waals (vdW-DF)®, and the range- 
separated hybrid (HSE)® xc-functionals have been carried out for 
comparison with the PBE based Go Wo calculations. In particular, 
we use the HSEO^^ variant of the HSE xc-functional. 

In the QP Go Wo approach, the contribution to the Kohn-Sham 
(KS) eigenvalues from the exchange and correlation (xc)-potential 
Vxc is replaced by the self energy E = /GW, where G is the Green’s 
function and W is the screening®! based on the KS wavefunc¬ 
tions.® The dielectric function is obtained from linear response 
time-dependent (TD) density functional theory (DFT) within the 
random phase approximation (RPA), including local field effects.® 
From Go Wo one obtains first-order QP corrections to the KS eigen¬ 
values, but retains the KS wavefunctions. Since our aim is to com¬ 
pare the computed interfacial level alignment with measured UPS 
spectra, it is most consistent to align the QP Go Wo levels with the 
VBM. 

We find Eyac, i-e., the effective potential far from the surface, 
from Go Wo is essentially the same as the Evac from DFT. In other 
words, the effective potential is unchanged by GqWq. To obtain 
a more accurate absolute QP level alignment relative to Evac, we 
employ a self-consistent QP GW approach.® In particular, by em¬ 
ploying the scQPGWl approach, we obtain both a QP PDOS com¬ 
parable to tha t from QP Go Wo and an improved alignment relative 
to Here, 25 %, 25 %, and 50 %, of the QP self energies 

are “mixed” with the DFT xc-potential over three self-consistent 
QP GW cycles,® respectively. If, instead, 100 % of the DFT xc- 
potential were replaced by QP self energy in a single self-consistent 
QP GW cycle, one would exactly obtain the QP Go Wo eigenvalues. 
However, this mixing is required to obtain a smooth convergence of 
both the QP wavefunctions and the absolute QP level alignment. To 
fully converge our self-consistent QP GW calculations (scQPGW), 
we perform a further eight cycles, with each introducing a further 
25 % of the QP self energy. 

The geometries have been fully relaxed using LDA,® PBE,® or 
vdW-DlSl xc-functionals, with all forces < 0.02 eV/A. HSE cal¬ 
culations are performed for the relaxed geometries obtained with 
PBE. We employ a plane-wave energy cutoff of 445 eV, an elec¬ 
tronic temperature ksT 0.2 eV with all energies extrapolated to 
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r ^ 0 K, and a PAW pseudopotential for Ti which includes the 3 s 
and 3 p^ semi-core levels. All calculations have been performed 
spin unpolarized. 

For the clean stoichiometric Ti02(l 10 ) surface® we have used a 
four layer slab and an orthorhombic 1x1 unit cell of 6.497 x 2.958 x 
40 A^i.e., 

/ V2a 0 0 \ 

0 c 0 , (1) 

\ 0 0 Vla + D ) 

where D 27 A is the vacuum thickness and a and c are the ex- 
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perimental lattice parameters for bulk rutile Ti 02 (a = 4.5941 A, 
c = 2.958 A).^ We have employed a F-centered 4 x 8 x 1 k-point 
mesh, and 320 bands = 9!/3 unoccupied bands per atom, i.e. includ¬ 
ing all levels up to 26 eV above the valence band maximum (VBM). 

For the clean reduced Ti 02 _ 1 / 4 ( 1 10) surface we have used a mon¬ 
oclinic 1x2 unit cell of 6.497 x 5.916 x 40 A^, i.e., 

/ \/2a c 0 \ 

0 2 c 0 , ( 2 ) 

\ 0 0 V2a + D ) 


to maximize the separation between the For the H 2 O covered 
surfaces, we have employed a four layer slab with adsorbates on 
both sides and an orthorhombic 1x2 unit cell of 6.497 x 5.916 x 

47 A^i.e., 

/ Via 0 0 \ 

0 2c 0 ’ (3) 

\ 0 0 V2a + D ) 


where D ^ 34 A. We employed a F centered 4x4x1 k-point mesh, 
with approximately We unoccupied bands per atom, i.e. including 
all levels up to 30 eV above the VBM, an energy cutoff of 80 
eV for the number of G-vectors, and a sampling of 80 frequency 
points for the dielectric function. The GqWq parameters are con¬ 
sistent with those previously used for descr ibing rutile Ti 02 bulk, 
Ti02(l 10) clean surface and interfaces.' —'— 'These parameters have 
been shown to provide accurate descriptions of bulk optical ab- 
sorpt ion sp ectra, and both clean surface and interfacial level align- 

ment.iaS 

To model H 2 O in the gas phase, we employed a unit cell with 
C 2 v symmetry and 16 A of vacuum in each direction. At the GqWq 
level, we used a smaller energy cutoff of 40 eV for the number of G- 
vectors, which has previously shown to provide an accurate descrip¬ 
tion of the optical absorption spectra for isolated molecules. 

To obtain DFT total energies and the relaxed structure of the 
clean reduced TiO 2 _i/ 8 ( 110 ) we have used a monoclinic 1x4 unit 


cell of 6.497 x 11.832 x 28 A^ i.e., 

/ ^/2a 2c 0 \ 

4c 0 , (4) 

\ 0 0 V2a + D ) 


where Z) ^ 15 A, and employed a F-centered 4x2x1 k-point mesh. 

In this study, we have performed PBE and subsequent single¬ 
point RPBE^ based DET calculations for the H 2 O adsorption en¬ 
ergies Eads on the stoichiometric and reduced surfaces. The RPBE 
xc-functional was especially developed for the prediction of adsorp¬ 
tion properties on metal surfaces.^ The H 2 O adsorption energy on 
the Ticus site of a stoichiometric TiO 2 ( 110 ) surface is given by 

g["H2O.T1O,(110)l-E[T,0,(110)] 

n 

where n is the number of adsorbed H 2 O functional units in the su¬ 
percell, and£[wH 2 O-rTiO 2 ( 110 )], £[TiO 2 ( 110 )], andE[H 20 ] are 
the total energies of the covered and clean stoichiometric surfaces 
and gas phase water molecule, respectively. Similarly, the H 2 O ad¬ 
sorption energy on the site of a reduced TiO 2 -jc( 110 ) surface 
is given by 


E'ads 


E[rM20 + Ti02-;,(110)] -£[Ti02-;,(l 10)] 
n 


-E[H20], (6) 


where Z^[wH 20 -r TiO 2 -jc( 110 )] and E[TiO 2 -jc( 110 )] are the total 
energies of the covered and clean reduced surfaces, respectively. 
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Figure 1. Intact H 2 O adsorbed with parallel (=^) interfacial hydrogen bonds 
on coordinately unsaturated Ti sites (H20@TicM5). (a) GqWq DOS for 1 ML 
of intact H 2 O covered (turquoise regions) or clean (gray region) stoichio¬ 
metric TiO2(110), their total DOS difference (dashed line), and the H 2 O 
PDOS. (b) Selected molecular orbitals at F and their energies (dotted lines). 
UPS difference spectra for H 2 O covered TiO2(110) (c) after 0.2 L expo- 
sure for T = 150, 160, 175, and 190 KPand (d l for r = 160 K after 0.05, 
0.1, 0.3, 0.7, and 1 L exposure.li^Peak position^^^^ are marked in brown, 
(e) H 2 O molecular orbitals, GqWq calculated eigenenergies marked in cyan, 
and experimental gas phase spectrum aligned with the Ibi level of (c).l^ 
Energies are relative to the VBM (^vbm)- Intensity references are provided 
for E > eybm when available. 


3. RESULTS AND DISCUSSION 

3.1. Intact H 2 O on the Stoichiometric Surface. In Figure]^ 
we disentangle adsorbate and substrate contributions to the spec¬ 
trum of intact H20@TicM,y, and compare the H 2 O PDOS to the the¬ 
oretical and experimental difference DOS. Specifically, we model a 
monolayer (ML) of H 2 O molecules with parallel (=^) interfaci al hv- 
drogen bonds aligned along the [001] direction (Figure [ijb)).^^^^ 
Note that IML of intact H 2 O is the most stable coverage and struc¬ 
ture on the stoichiometric rutile TiO 2 ( 110 ) surface.^ 

The theoretical difference DOS is the difference between the to¬ 
tal DOS of the H 2 O covered (H20@TicM,y) and clean stoichiomet¬ 
ric (TiO2(110)) surfaces, as shown schematically in Figure [^a). 
Turquoise areas in the ll20@Ticus and difference DOS indicate re¬ 
gions of greater density for the H 2 O covered versus clean stoichio¬ 
metric surface. The gray area indicates the DOS energy range for 
the clean stoichiometric TiO2(110) surface. Figure [^c) and (d) 
show two sets of UPS difference spectra obtained either by raising 
the temperature (from 150 K to 190 K) for a consistent ex^sure to 
H 2 O (0.2 L) for an annealed TiO2(110) surfacel^(Figure^c)), or 
by increasing the H 2 O dose (from 0.01 L to 1 at low temperature 
(160 K) for a nearly perfect surface!^ (Figure [^). The experimen¬ 
tal spectra have been referenced to the VBM, which is positioned 
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3.2 eV below the experimental Fermi level.^ 

Comparing the difference DOS to the H 2 O PDOS, we find the 
peaks lying outside the TiO2(110) DOS energy range are clearly 
attributable to H 2 O levels. As shown in Figure [^b), these levels 
are related to the lb 2 and 3ai H 2 O orbitals shown in Figure [^. 
This is not the case within the TiO2(110) DOS region, where the 
adsorbate levels are broadened by hybridization with the surface. 
This hybridization with the surface has been severely underesti¬ 
mated by previous cluster-based MP2 calculations.^ Within the 
TiO2(110) DOS region, the peaks in the H 2 O PDOS have corre¬ 
sponding peaks in the difference DOS, although the relative peak 
intensities differ substantially between the two methods. More im¬ 
portantly, the difference DOS has dips centered at -4.1, -2.4, and 
-1.1 eV, where there are adsorbate levels in the PDOS, and a peak 
at -0.4 eV, where there are no adsorbate levels in the PDOS. The 
dips at -4.1 and -1.1 eV correspond to the O 2po- and O 2p;r ^aks 
in the TiO2(110) DOS, ^ respectively, as marked in Figurej^a). 
These peaks split due to mixing with the 3ai and Ibi H 2 O orbitals. 
This splitting is the origin of the observed dips in the difference 
DOS, which are also seen experimentally in Figure [^c) and d. 

The peak at -9.4 eV in the H 2 O PDOS, which has lb 2 molec¬ 
ular character, agrees semi-quantitatively with the most strongly 
bound experimental peaks at -9.8 eV (Figure [^c)) or -10.0 eV 
(Figure [^). The peak at -6.3 eV in the H 2 O PDOS, which has 
intermolecular 3ai bonding character, agrees semi-quantitatively 
with the experimental peaks at -6.4 eV (Figure [^c)) or -7.1 eV 
(Figure [^). Note that the theoretical average deviation is within 
that amongst the experiments. This may reflect differences in sam¬ 
ple preparation, which result in a variety of different H 2 O config¬ 
urations, i.e., H 2 O coverages, O^^^ concentrations, and mixtures 
of intact and dissociated H 2 O. As we will show in Section [34] by 
considering a variety of H 2 O structures a more complete descrip¬ 
tion of the experiment is obtained. Altogether, this agreement for 
the -9.4 and -6.3 eV PDOS peaks lends confidence to our results 
for regions where the experimental results are unclear. 

The assignment of the peaks located within the TiO2(110) DOS 
is much more complicated. The assumption that the highest peak 
in t he exp erimental spectra originates solely from the H 2 O Ibi 
level^^^^is an oversimplification. In fact, both the 3ai and Ibi 
molecular levels contribute within this region (Figure [^b)). While 
the levels with intermolecular 3ai bonding character give rise to a 
distinct peak below the TiO2(110) DOS region, those with inter¬ 
molecular 3ai antibonding character are pushed to higher energies 
and mixed with the Ibi molecular levels (Figure [^b)). The lat¬ 
ter is due to symmetry breaking at the interface. Consequently, the 
H 2 O PDOS is broadened into several peaks between -5 and -1 eV. 
These levels have interfacial (3ai/lbi- O 2po-/2p;r) bonding and 
antibonding character (not visible at the isosurface value used). 

3.2. Dissociated H 2 O on Reduced Surfaces. To see how 
dissociation of H 20 @ 0 ^^‘^ affects the spectrum, we now consider 
ViML of H 2 O dissociated on a reduced TiO2_i/4(110) surface (Fig¬ 
ure 0. Here, we have used TiO2_i/4(110) to denote a surface con¬ 
sisting of ViML of defects. This structure corresponds to the 
staggered surface species, shown in Figure [^b). 

The theoretical difference DOS is the difference between the to¬ 
tal DOS of the H 2 O covered (H20@0^^^) and the clean reduced 
(TiO2_i/4(110)) surfaces, shown schematically in Figure l^a). 
Turquoise areas in the H20@0^^ and difference DOS indicate re¬ 
gions of greater density for the H 2 O covered versus clean reduced 
surface. The gray area indicates the DOS energy range for the clean 
reduced TiO2_i/4(110) surface. The defects give rise to occu¬ 
pied levels with Ti 3d character that are just below the conduction 
band minimum and outside the energy range shown.^^^^^^Note that 
the H 2 O PDOS includes half the O atoms and all the H atoms that 
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Figure 2. H 2 O dissociated on bridging O vacancies (H 20 @ 0 ^^^). (a) 
GqWq DOS for V 2 ML of dissociated H 2 O covered (turquoise regions) or 
clean (gray region) defective TiO 2 _i/ 4 ( 110 ) with YiML of their total 
DOS difference (dashed line), and the H 2 O PDOS. (b) Selected molecular 
orbitals and their energies (dotted lines). UPS difference spectra for H 9 O on 
reduced TiO2-x(110) (c) after 0.2 L exposure for T = 260 and 300 Kpl[d) 
for T = 300 K after between 0.01 and 100 L exposure, EH and (e) for T = 
120 K after 0.14, 0.3, 0.4, 0.5, and 0.7 L exposure .EH Peak positionsE^EHHl 
are marked in brown. Energies are relative to the VBM (sybm)- Intensity 
references are provided for s > ^vbm when available. 


make up the species. In this way the PDOS is provided in 
terms of H 2 O formula units. 

The peak in the difference DOS and PDOS at -7.0 eV has 0/,^H 
cr character, as shown in Figure |^b). Note that the peak intensity 
in the PDOS is about half that in the difference DOS, as the PDOS 
includes half the O^r atoms. This peak’s position agrees semi- 
quantitatively with the experimental peaks at -7.1 (Figure [^c)), 
-7.6 (Figurej^), or -7.2 eV (Figurej^). The PDOS has a broader 
feature between -4 and -1 eV, due to hybridization with the sur¬ 
face. This feature is associated with contributions coming from the 
bonding and antibonding combinations of two distinct p orbitals of 
the OhrU species (Figure |^b)): one per pen dicular to the 0/,^H cr 
bonds (the so-called OH n level of NaOH^; the other in the plane 
of the OhrU cr bonds. The lowest of these peaks at -3.9 eV cor¬ 
responds to the bonding combination of the 0/,^H n levels. This 
peak’s position agrees semi-quantitatively with the consistently ob¬ 
served experimental peaks at -4.8, -4.4, and -4.5 eV in Figures 
Scxiiji, and|^, respectively. However, the antibonding 0/,^H n 
levels are shifted to much higher energies (-1.2 eV), as shown in 
Figurej^b). 

Much of the theoretical difference DOS’s structure is attributable 
to the defect healing of as seen from the difference DOS be¬ 
tween TiO2(110) and TiO2_i/4(110) in Figurel^ This suggests that 
the observed features in the experimental difference spectra over- 
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Figure 3. difference DOS between (red) stoichiometric Ti 02 (l 10) and 
(black) reduced Ti 02 _ 1 / 4 ( 1 10) with ViML of defects. Red areas indicate 
defect healing of i.e., regions of greater density for the stoichiometric 
versus reduced surfaces, shown in Figures|^and[^ respectively. 




Figure 4. (a) IML intact H 2 O adsorbed with parallel (=^) interfacial hy¬ 
drogen bonds (black dashed lines) on coordinately unsaturated Ti sites 
(H 20 @TicM 5 ). Total (maroon) and H 2 O projected (blue) DOS computed 
with (b,d,f) DFT and (c,e,g) GqWq using the (b,c) local density anproxima- 
tion (LDA)E^ (d,e) generalized gradient approximation (PBE)I^ and (f,g) 
long-ranged van der Waals interactions (vdW-DF)l^ for the xc-functional. 
The calculated H 2 O PDOS are compared with the UPS spectrum at 150 K 
after 0.2 L exposureEH (black). Energies are relative to the valence band 
maximum, ^vbm- 


Table 1. Height of H 2 O Above TiO2(110) for IML Intact H20@Tic„s 
Measured with SXPS and Calculated with LDA, PBE, or vdW-DF XC- 
Functionals. 


Method 

d[H20-Ticus] (A) 

sxpsi^^ 

2.210 

LDA 

2.180 

PBE 

2.367 

vdW-DF 

2.434 


lapping with the reduced surface’s DOS are simply levels rein¬ 
troduced by dissociated H 20 @ 0 ^^^. In particular, the peak which 
is usually attributed to n levels is actually composed of 
surface levels unrelated to the presence of H atoms. 

3.3. XC-Functional and Methodology Dependence of H 2 O 
Spectra for Stoichiometric and Reduced Surfaces. To assess 
the robustness of the calculated QP H 2 O PDOS, we consider its 
dependence on the xc-functional and methodology. Specifically, we 
compare the H 2 O PDOS from DFT, scQPGWl, and GqWo for IML 
intact ll 20 @Ticus with parallel (^) and antiparallel (?±) interfacial 
hydrogen bonds and ViML dissociated H 20 @ 0 ^^^ in Figures 
an respectively. 

We find the observed structure of the GqWq H 2 O PDOS is in¬ 
dependent of whether the local density approximation (LDA),^ 



Figure 5. (a) IML intact H 2 O adsorbed with antiparallel (?^) interfacial 
hydrogen bonds on coordinately unsaturated Ti sites (H 20 @TicM 5 ). To¬ 
tal (maroon) and H 2 O projected (blue) DOS computed with (b,c) DFT, 
(d) scQPGlTl, and (e.f) GqWq using the (b,d,e) generalized gradient ap¬ 
proximation (PBE)^ and (c,f) range-separated hybrid (HSE06)E1 for the 
xc-functional. The calculated H 2 O PDOS are compared with the UPS spec¬ 
trum at 150 K after 0.2 L exposure^H (black). Energies are relative to the 
valence band maximum, ^vbm- 



DFT 



Figure 6. (a) ViML H 2 O dissociated on bridging O vacancies (H 20 @ 0 ^^^) 
of defective Ti 02 - 1 / 4 ( 1 10) with / 2 ML of O^^^. Total (maroon) and H 2 O pro¬ 
jected (blue) DOS computed with (b,c) DFT and (d,e) Go Wq using the (b,d) 
generalized grad ient approximation (PBE)ISland (c,e) range-separated hy¬ 
brid (HSEOb/^for the xc-functional. Filling denotes occupation. The cal¬ 
culated H 2 O PDOS are compared with the UPS spectrum at 300 K after 0.2 
L exposureEl (black). Energies are relative to the valence band maximum, 

^VBM- 


generalized gradient approximation (PBE),I^ long-ranged van 
der Waals interactions (vdW-DF)^, or a range-separated hybrid 
(HSE06)E3 are employed for the xc-functional. This is consistent 
with the previously reported similarities between PBE and HSE 
based Go Wo PDOS for CH 3 OH on TiO 2 ( 110 )PThis is despite 
the greater differences observed amongst the DET H 2 O PDOS, 
which all differ qualitatively from the experiments. Furthermore, 
the GqWq H 2 O PDOS is robust to the resulting changes in the 
H 2 O height above the surface, i.e., the distance between H 2 O and 
Ticus d[}l20 -Ticus], shown in Table[^ Furthermore, Figure[^d,e) 
shows that scQPGWl provides a similar H 2 O PDOS level aligne- 
ment to GqWq. This is consistent wit h what was previously reported 

for the CH 3 OH-TiO 2 ( 110 ) interface.^Sl 

We clearly see that the differences between the DFT and GqWq 
PDOS, i.e., the QP energy shifts, are far from simply being rigid. 
For instance, we find for PBE that the QP energy shifts for the 
levels that contribute to the highest-energy PDOS peak are 

almost negligible (c/ Eigures |^d,e) and[^b,e)) . As a result, the 
QP GqWq is only ~ 0.1 eV lower compared to DPT. On the 

other hand, we find significant QP shifts to stronger binding for the 
levels that contribute to the most strongly bound PDOS peak with 
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Ib2 cr molecular character. For example, with PBE the QP GqIFo 
lowest energy peak is shifted by ~ - 1.7 eV compared to DFT (cf. 
Figures|^d,e) and[^b,e)). 

As previously shown for the CH30H-Ti02(l 10 ) interface, these 
differences in the shifts of the peaks are directly related to differ¬ 
ences in the spatial distribution of the wave functions for the levels 
contributing to the peaks.*—" — 'This is because the QP GqWq cor¬ 
rections to the DFT eigenenergies for interfaces are directly corre¬ 
lated with the spacial distribution of the wave functions.The 
negligible shift of the DFT highest-energy PDOS peak (Figures 
(b,d,f) and|^b,c)) is due to its strong hybridization with the sur¬ 
face, i.e., weight on TiO2(110), for the levels contributing to this 
peak.^^^’^On the other hand, the levels that contribute to the most 
strongly bound PDOS peak have little weight on TiO2(110), and 
have cr character. Both their localized H2O character as well as 
their cr nature explain why these levels have large QP energy shifts 
to stronger binding. 

Oxygen defective and hydroxylated (/i-)Ti02 surfaces have oc¬ 
cupied 3 d levels which are associated with reduced Ti^"^ atoms.^ 
One such example is the ViML dissociated H20@(^^ on reduced 
TiO2-i/4(110) with ViML of O^^^ shown in Figure [^a). The spa¬ 
cial distribution of the 3 d density for O defective surfaces has been 
charact erized by low temperature scanning tunneling microscopy 
(STM)STM measurements find at 77 K the 3 d density is ho¬ 
mogeneously distributed along the [ 001 ] direction,^ while at ~ 5 K 
the 3 d density exhibits an asymmetric localized character.^ 

A localized description of the Ti^'*' occupied 3 d levels is not ob¬ 
tained from DFT with standard xc-functionals. For example, the 
occupied 3 d levels obtained with PBE are highly delocalized, as 
clearly shown in Eigurej^b). This is due to self-interaction errors 
which are inherent in such xc-functionals. If one performs spin- 
polarized DET calculations with a hybrid xc-functional on such 
systems, one obtains localized Ti^"^ 3 d^ levels between 0.7 and 
1.6 eV below the C BM a long with a structural deformation of 
the TiO2(110) surf ace.^^^® However, spin-paired calculations with 
HSE 06 on the PBE relaxed geometry only yield an occupied shoul¬ 
der at the CBM (Eigure |^c)). At the QP GqWq level based on 
PBE, this shoulder evolves into a distinct peak about 0.6 eV below 
the Eermi level, ep- This effect is even more pronouced when the 
Go Wo calculation is based on HSE 06 (cf. Eigure [^d,e)), which 
yields peaks at 0.6 and 0.9 eV below ep- As compared to Go Wo 
PBE, Go Wo HSE 06 shifts the unoccupied 3 d levels further up in 
energy revealing the double peak structure. These energies are in 
very good agreement with the peak at 0.8 eV below ep in the UPS 
spectra of H20@0^^^ of Eigurej^d). This peak is not shown in 
Figure |^d) as it is slightly above 2 eV with respect to VBM.^ 
However, note that Go Wo overestimates by about 1 eV the VBM 

1101 

position relative to &p as compared with UPS experiments.' —' 

This result is completely independent of the wavefunction’s spa¬ 
cial distribution, i.e., localization, as the Go Wo calculations are 
based on the KS wavefunctions. This is different from previ¬ 
ous findings, which showed DPT with either PBE or hybrid xc- 
functionals is only giving distinct peaks for the occupied 3 d levels 
provided the relax ed spin-polarized distorted structure is used in the 
calculations 

While for Go Wo based on PBE and HSE 06 one sees noticeable 
differences in the description of the 3 d occupied levels, the QP 
H2O PDOS and its alignment relative to the VBM are unchanged. 
Although localization of the Ti^"^ occupied levels and associated 
structural deformations are absent from our approach, such features 
should not significantly alter the QP H2O PDOS. This is because 
the Ti^"^ levels are too far above the VBM (~ 2 eV^^ to hybridize 
with the H2O. Moreover, as we will show in Section |T 4 | the QP 
H2O PDOS is rather robust to local deformations of the surface 


structure, e.g., due to changes in coverage. 

3.4. Coverage and Dissociation Dependence of H 2 O Spec¬ 
tra for Stoichiometric and Reduced Surfaces. As different 
experimental conditions and surface preparations have been em¬ 
ployed, there are expected to be different H2O structures on the 
surface. To evaluate how strongly the DOS depends on the adsorp¬ 
tion geometry, we now consider a variety of coverages of intact and 
dissociated H2O on rutile stoichiometric TiO2(110) (Eigure^ and 
reduced Ti02_ 1 /4 ( 110 ) (Eigurej^ and Ti02_ 1 /g ( 110 ) (Eigurej^ with 
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Dissociation 


Figure 7. Schematics of H 2 O adsorbed intact (I) or dissociated (D) on co- 
ordinately unsaturated Ti sites (Tic^^) of stoichiometric TiO2(110). Higher 
coverages are obtained by the addition of second-layer H 2 O. Coverage is 
the number of H 2 O formula units per (110) 1x1 unit area of the clean stoi¬ 
chiometric surface. Dissociation is the fraction of H 2 O molecules which are 
dissociated, i.e., one minus the ratio of intact H 2 O molecules to H 2 O for¬ 
mula units. Colored frames encompass regions of common fractional disso¬ 
ciation. Charge transfer of about -0.4e accompanying deprotonationl^ of 
intact H 2 O adsorbed at Tic^^ is represented by arrows, while intermolecular 
(gray) and interfacial (black) hydrogen bonds are denoted by dotted lines. 


ViML and V4ML of Ojff defects, respectively. The relative impor¬ 
tance of these geometries is illustrated in Eigurep^a) andp^b) by 
the average absorption energyper H2O molecule on the stoi¬ 
chiometric or reduced surfaces®with either PBE^or rpbeP xc- 
functionals. In so doing, the contribution of different structures to 
the measured spectra can be disentangled. Note that an intact V2ML 
of H20@0^^^ (Eigurej^b)) is probably only a transient locally 
stable state of the reduced H20-Ti02-1/4(110) interface,!^ which 
may easily evolve into the ~ 0.7 eV more stable dissociated V2ML 
H20@0^^^ (Eigure|^c)). Eor this reason, we only consider disso¬ 
ciated H20@0^^‘^ structures in Figure p^. , 

By co mparin g to lower coverage H9O structures 
to lML™ 22 | 621 in Figure and ‘^ML^in Figure § to m 

Figure [^, we can disentangle the effect of interaction between the 
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Figure 8. Schematics of reduced TiO2_i/4(110) with ViML of bridging O vacancies (a) clean, covered with YiML (b) intact and (c) dissociated 

H20@0^^^, and with an additional (d) ViML or (e) IML of intact H 2 O adsorbed on coordinately unsaturated Ti sites (Ticus)- Coverage is the number of H 2 O 
formula units per (110) 1 x 1 unit area of the clean reduced surface. Dissociation is the fraction of H 2 O molecules which are dissociated, i.e., one minus the ratio 
of intact H 2 O molecules to H 2 O formula units. Charge transfer of about -OAe accompanying deprotonationl^is represented by arrows, while intermolecular 
(gray) and interfacial (black) hydrogen bonds are denoted by dotted lines. 
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Figure 9. Schematics of H 2 O adsorbed dissociated (D) on V 4 ML of bridg¬ 
ing O vacancies (O^^^) on reduced TiO2-x(110) (x = Vs). Higher coverages 
are obtained by the addition of H20@TicM5. Coverage is the number of 
H 2 O formula units per (110) 1x1 unit area of the clean stoichiometric or 
reduced surface. Dissociation is the fraction of H 2 O molecules which are 
dissociated, i.e., one minus the ratio of intact H 2 O molecules to H 2 O for¬ 
mula units. Colored frames encompass regions of common fractional dis¬ 
sociation. Charge transfer of about -OAe accompanying deprotonationl^ 
of intact H 2 O adsorbed at Tic^^ or is represented by arrows, while in¬ 
termolecular (gray) and interfacial (black) hydrogen bonds are denoted by 
dotted lines. 


H2O molecules on the spectra. Further, these structures allow us to 
probe the isolated molecule limit. 

As shown in Figure at lower coverages the overall width of 
the spectra is reduced with fewer distinct peaks. When the coverage 
is increased to include intermolecular interactions between adjacent 
species, the molecular levels hybridize into bonding and antibond¬ 
ing intermolecular levels. This produces additional peaks above 
and below those present at low coverage. As a result, the peak 
with intermolecular bonding 3 ai character at - 6.3 eV for IML of 


ll 20 @Ticus is absent for a V2ML coverage. This reinforces the as¬ 
signment of the experimental spectra shown in Figure[^to an intact 
IML ll 20 @Ticus geometry with interacting molecules. 

To see how the spectra for dissociation of ll 20 @Ticus compare 
to H20@0^^^, we have considered the half-dissociated (ViD) and 
fully dissociated (D) H2O structures shown in FigureAs shown 
in Figure p^c), the peak at - 7.0 eV with cr character for 

H20@0^^‘^ splits into two peaks for dissociated ll 20 @Ticus- The 
lower energy peak has both and cr character, while 

the higher energy peak is mostly in character. Furthermore, 

we find a similar couple of peaks for y4ML mixtures of dissociated 
}l 20 @Ticus and H20@0^^‘^ shown in Figure [T^. This means one 
may recognize dissociated ll20@Ticus by both the presence of two 
peaks at about - 7.0 and - 6.3 eV, and the absence of the low-energy 
peak with lb2 character for intact ll20@Ticus’ 

The absence of a peak at about - 6.3 eV in the experimental 
spectra shown in Figure |^c) reinforces its attribution to dissoci¬ 
ated H20@0^^^ rather than dissociated H20@TicM5. This is further 
supported by the calculated H2O absorption energies (Figure p^a) 
andp^b)). These are generally weaker for dissociated H 90 @TWv, 
and stronger for H20@0^^^, as in previous calculations.^ 

To check whether changes in the absorption geometry of H2O af¬ 
fect the spectra for the same coverage, we compare IML of H2O {I, 
V2D, D} adsorbed with either parallel (=^) or antiparallel (?±) inter¬ 
facial hydrogen bonds^ (black dashed lines in Figure]^. Overall, 
the two sets of spectra are consistent, and demonstrate the general 
robustness of the DOS to minor changes in the water absorption 
geometry. However, as the H2O molecules are no longer equiva¬ 
lent when the interfacial hydrogen bonds are antiparallel, there is 
a greater splitting between bonding and antibonding contributions 
for the peaks with lb2 and 3 ai molecular character. In particular, 
for intact H2O, the lowest energy peak with molecular lb2 char¬ 
acter splits with a separate peak at - 9.6 eV, which is closer to the 
peaks at - 9 . 8 ^(Figure[^c)) and - 10.0 eV® (Figure[^) observed 
experimentally. 

To see how increasing the H2O coverage impacts the spectra, we 
compare monolayer (WML or IML) to multilayer (IWML) H2O 
{I, WD, (Figure [^, and consider the effect of additional 

H20@Tic„i to ‘/4ML (Figure!^ and ViML (Figure]^ Yi 20 @OlfW^ 
In this way we can can see how robust the observed features in 
the individual spectra for isolated species are to screening by H2O 
layers,® and probe the liquid water limit.^ 

When a second layer of H2O is added to the low coverage intact 
y2ML H20@TicM.y structure, the levels with H2O lb2 character are 
unchanged, while the levels with 3 ai and Ibi second layer character 
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Figure 10. Structure and coverage dependence of (a,b) adsorption energy Eads and (c,d) GqWq PDOS for H 2 O adsorbed intact (I) or dissociated (D) on (a,c) 
coordinately unsaturated Ti sites (Ticus) of stoichiometric TiO2(110) (Figurel^ and (b,d) bridging O vacancies of reduced TiO2-x(110), with x = Vs 

(thin lines, Figurej^ or y4 (thick lines, Figurejsj. (a,b) Eads calculated with PBE (O) and RPBE (O) xc-functionals for (white) low (y4 and y 2 ML), (turquoise) 
medium (y4 and IML), and (blue) high (IV4 and IY2ML) coverage. UPS difference spectra at (c) 150 K and (d) 300 K after 0.2 L exposure are from Ref. ED 
(c,d) Energies are relative to the VBM (evBivi)- Gray regions denote the clean surface DOS. Red dashed lines denote the highest PDOS peaks for IML 

H20@TW and '/ 2 ML H20@0^“. 


are more localized and weakly hybridized with the surface. These 
levels are seen as the two most intense peaks at -4.3 and -2.2 eV 
(Figure p^c)). The former coincides with the peak at -4.2 eV ob¬ 
served experimentally at low temperatures (Figure[2c)), suggesting 
multilayer H2O structures may be present under these experimental 
conditions. The intermolecular H bonding between the layers delo¬ 
calizes the molecular levels of the first layer. This is seen from the 
peak at -6.1 eV with antibonding 3ai character on the first layer. 
We saw the same behavior when increasing the first layer’s cover¬ 
age from y2ML to 1 ML. This is further confirmation that the peak 
observed experimentally at -6.4 eV has intermolecular character. 

When a second y2 layer of H2O is added to the IML H20@TicM5 
{UsD, ysD} structures,^ a denser network of intermolecular and 
interfacial hydrogen bonds is formed, as shown in Figure]^ This 
causes a stronger hybridization between the OH and H2O cr levels. 
For the UsD structure, this results in the four distinct cr peaks shown 
in Figure p^c). On the one hand, the peaks at -9.1 and -6.2 eV 
have predominantly intact H2O and Ocus^ character, as was the 
case for IML of y2D H20@TicM.y. On the other hand, the peaks at 
-7.9 and -7.4 eV are most related to the second layer. In effect, the 
H2O cr level of the second-layer H2O, which is fully saturated with 
four hydrogen bonds, is upshifted by more than an eV. 

This is not the case for the ViD structure (Figure [^, where the 
peak at -9.1 eV instead has mostly intact second-layer H2O lb2 
character. As was the case for intact lYiML H20@TicM.y, the ad¬ 
dition of a second 1/2 layer of H2O induces a stronger hybridiza¬ 
tion of the Ohrii levels, and introduces an additional intense peak 
at -4.4 eV (Figure p^c)). This again suggests the experimentally 
observed peak at -4.2 eV (Figure ^c)) may be due to multilayer 
H2O. 

Overall, we find the addition of second-layer H2O affects the re¬ 
sulting spectrum qualitatively. We find both additional features and 
a redistribution of those due to the first H2O layer. When we in¬ 
stead add H20@TicM5 to the V4ML and V2ML H20@0^^‘^ structures 
(Figures]^ and[^ we find the resulting spectrum is the sum of the 
separate spectra to within 0.2 eV (Figure [^. For example, the 
U/2ML UsD spectrum (Figure p^) for IML of intact H2O added to 
V2ML H20@0^^‘^ (Figure]^ is basically the sum of the IML in¬ 


tact H20@Tici,, (Figurel^a)) and ViML H20@0^^^ (Figure |2];a)) 
PDOS spectra downshifted by 0.2 eV. This explains the ease with 
which the experimental single-layer H2O spectra may be analyzed 
for levels outside the surface DOS region. 

3.5. Alignment of the Highest H 2 O Occupied Levels. So 
far, we have concentrated our analysis on the lower energy peaks 
observed in the experimental spectra. This was done to demonstrate 
the robustness of the calculated QP DOS. Having established this, 
we now focus on the adsorbate levels near the VBM, which play 
an important role in photooxidation processes. In this respect, the 
highest H2O occupied levels’ alignment for IML intact and dis¬ 
sociated H20@TicM5, and V2ML dissociated H20@0^^^ is of ut¬ 
most importance. The former structure corresponds to the reac¬ 
tant species on stoichiometric surfaces,^ which undergoes photo¬ 
irradiation. The latter structures act as ho le traps and are thus the 
main oxidizing agents on TiO2(110).^^^^ 

We have shown that the experimental peak at -4.2 eV^is not, 
in fact, the highest energy peak of H20@TicM.y. We instead find 
the highest-energy PDOS peak, for IML intact H20@TicM.y 

at -1.1 eV relative to the VBM (Figure p^c)). This is 0.6 eV 
closer to the VBM than the —1.7 eV estimate!^ deduced from 
the onsets of the UPS difference spectra in Ref. [18] Moreover, as 
IML H20@TicM.y dissociates, moves up to -1.0 eV (ViD) 

and -0.9 eV (D) (Figure [T^)). This is again significantly higher 
than the —1.8 eV estimate^ based on UPS difference spectra for 
the TiO2(100) surface from Ref. |7^ As was the case for CH3OH 
on TiO2(110),^this raising of Sp^^^ can be related to the charge 
transfer of -0.4^ that accompanies deprotonation (arrows in Fig¬ 
ure [^. We find for the IML intact structure on Ti02(l 10) Sp^J^^ is 
0.2 eV closer to the VBM for H 2 O than for CHaOH.ESESwhile for 
the IML '/ 2 D structures is the same. However, the highest 

PDOS peak is both less intense and broader for H2O compared to 
CH3OH, due to the stronger hy bridi zation with the surface. This 
is why, as discussed in Section |3.l[ the QP GqWq is only 

~ 0.1 eV lower compared to DFT^^^^^ (Figure]^. After adding 
second-layer H2O, is unchanged with weight mostly remain¬ 

ing on the first layer. 

We find for VML dissociated H20@0^“^ r; -1.1 eV rel- 
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ative to the VBM (Figure p^), the same as for intact ll 20 @Ticus- 
This is much higher than the previous estimate of ~ - 3.7 eV^for 
Ohrii based on the UPS difference spectra in Ref. Our cor¬ 
rected value agrees with the recently demonstrated photo- 

catalytic inmortance of sites as the main oxidizing species on 

TiO2(110).l^ 

Based on for IML intact H20@TicM5, vertical excitations 

from the highest H2O occupied levels to the TiO2(110) conduction 
band require photon energies that exceed the electronic band gap 
for bulk rutile Ti02 ( 3.3 ± 0.5 eV^ by > 1 eV. However, the hole 
generated by such supra-band gap excitations should be mostly 
located on TiO2(110) O 2 p;r rather than H2O O 2 p levels. This is 
because the H2O highest levels are hybridized with TiO2(110) and 
are predominantly TiO2(110) in character. 

The fact that the highest H2O levels are ~ 1 eV below the VBM 
does not necessarily mean that they cannot be photooxidized by 
holes photogenerated within the TiO2(110) valence band. A recent 
DFT study with HSE 06 found trapped holes at surface O sites, i.e., 
three-fold coordinated Osfaid, are shared with nearyby HO-TicM^ 
groups.^ Moreover, it has been suggested that H 9 O can only be 
photooxidized, i.e., trap a hole, upon deprotonation.' —'—' In other 
words, hole transfer to the HO-TicM^ site should be mediated by 
the deprotonation of intact H20@TicM5 to the nearest O^r site. Al¬ 
together, this suggests that ll20@Ticus photooxidation should be 
initiated by band-to-band and supra-band photo-excitations, which 
result in the generation of holes within the Ti02(l 10) valence band. 
These Ti02(l 10 ) free holes may then be trapped at sites, and 
partially transferred to nearby HO-TicM^ upon H2O deprotonation. 

3.6. Vacuum Level Alignment. So far, we have considered the 
level alignment of the interfacial levels relative to the VBM of the 
substrate. This allows a direct comparison of the occupied PDOS 
with the measured UPS spectra. However, to assess the photoelec- 
trocatalytic activity of the interface, one needs the absolute level 
alignment relative to the vacuum level Evac- 

In Figure[TT]we show the level alignment for gas phase H2O and 
IML intact H20@TicM5 relative to E^ac from DFT, scQPGIVl, and 
Go Wo based on PBE and HSE xc-functionals. These are com pared 
to the measured CBM for the liquid H20-Ti02(l 10 ) interfacel^® 
and the measured and coupled-cluster (CCSD(T)) gas phase H2O 
ionization potential.^ 

Our calculated IP values for H2O in gas ph ase y e consistent 
with those reported previously in the literature } Although 
the relative energies of the Ibi, 3 ai, and lb2 H2O levels are con¬ 
sistent over all five levels of theory, the levels are rigidly down¬ 
shifted. We observe a clear ordering in increasing IP of PBE DET 
( 7.2 eV) < HSE DET « PBE scQPGWl < PBE Go Wo < HSE Go Wo 
< PBE scQPGW ( 12.8 eV) < Hartree Fock (HE 13.9 eV®). 

To understand the origin of this ordering, we have probed the de¬ 
pendence of the IP on the fraction of Hartree-Eock exact exchange 
included in the range-separated HSE xc-functional via the param¬ 
eter a in Eigure On the one hand, for DET, we find a strong 
linear dependence of IP on a, i.e., IP /Ppbe + {lPa=\ “^^pbe)q^ 
7.2 -r 5 . 9 Qf, with a ~ 0.9 providing a quantitative agrement with 
experiment and CCD(T) calculations. Overall, this linear depen¬ 
dence is not surprising, as a may be interpreted as the amount 
of e lectro n-electron screening, i.e., the inverse dielectric constant 
,1221211 In other words, the fraction of exact exchange a included, 
determines the amount of screening, , incorporated within the 
xc-functional. The quantitative agreement of the IP for a ~ 0.9 is 
because small molecules, e.g., H2O, are weakly screened in the gas 
phase (^oo ~ 1). 

On the other hand, for GqWq, the calculated IP has a much 
weaker dependence on cr, i.e., the starting xc-functional, with IP 
/Ptt=i -A/P( 10 “”' - 1 ) Ri 13 . 4 - 1.2 X 10 ”“. Further, the GqWo and 
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Figure 11. Absolute level alignment for IML intact H2O adsorbed with 
antiparallel (?^) interfacial hydrogen bonds on coordinately unsaturated Ti 
sites (H20@TicM5). Total (maroon) and H2O projected (blue) DOS com¬ 
puted with DFT, scQPGlTl, and GqWq using the generalized gradient ap¬ 
proximation (PBE) and hybrid (HSE) xc-functionals. Energies are relative 
to the vacuum level Evac- The measured ^cbm from Ref. | 4 ^ (thick gray 
line), measured and coupled-cluster (CCSD(T)) H2O gas phase ionization 
potentials IP from Ref [|o|(thin gray line), and for each level of theory the 
calculated gas phase Ibi, 3 ai, and lb2 H2O levels (marked in cyan) are 
provided. 

Dielectric Constant 
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Figure 12. Ionization potential IP versus exact exchange fraction a in¬ 
cluded in the HSE xc-functional and equivalent dielectric constant £00 
from DFT (filled symbols), GqWq (open symbols), scQPGWl (green filled 
symbols), and scQPGW (brown filled symbols) for H2O in gas phase (cir¬ 
cles), a cleanly (diamonds) and a IML of intact H20@TicM5 (squares) on 
the stoichiometric Ti02(l 10 ) surface, and a ViML of dissociated H20@0^^^ 
(triangles) on the defective TiO2-i/4(110) surface with ViML of The 
mea sured IP for H2O in gas phase,the stoichiometric TiO2(110) sur- 
f^ce,lllllllillthe 6 - 9 % hydroxylated /z-Ti02(l 10 ) surface, ^^and the liquid 
H2O-TiO2(110) interface® are shown in gray. The self-consistent QP GW 
IP for H2O in gas phase is indicated by the horizontal dashed line. The 
experimental dielectric constant of bulk Ti02, « 7 . 6 ,® averaged over 

the ( 110 ) surface is marked in red. A linear fit to the DFT IP (blue), and an 
exponential fit to the GqWq IP (cyan) for H2O in gas phase are provided for 
comparison. 


DFT IP coincide when a ^ 1 . For GqWq based on PBE (a = 0 ), 
the IP already agrees semi-quantitatively with experiment, with 
full quantitative agreement obtained for GqWo based on HSE 06 
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(a = 0 . 25 ). This is because the RPA £oo ~ 1 , independently of 
a. Essentially, the calculated GqWq IPs would also be obtained 
from DFT using an HSE xc-functional with 0.84 < a < 1 . 0 , i.e., 
1 < ^oo < 1 . 2 . Overall, this implies GqWq is a predictive method for 
the IP of small molecules. However, the scQPGlT technique has 
the added advant age of being completely independent of the start¬ 
ing xc-functionalwhile providing a nearly quantitative IP. 

For the H2O-TiO2(110) interface, e.g., IML intact H20@TicM^, 
the highest energy H2O PDOS peak, , is pinned ~ 1 eV below 
the VBM across PBE DFT, HSE DFT, PBE scGWl, PBE GqWq, 
and HSE GqWq. For this reason, the IP of the H2O interfacial levels 
is controlled by the alignment of the VBM with respect to the vac¬ 
uum. This means we only need to consider the absolute VBM level 
alignment of the interface, i.e., the interface’s IP = -evBM + Evac, 
as a descriptor of photoelectrocatalytic activity. 

In Figure[^we see that the IP of the interface follows a different 
ordering across the methodologies from that of gas phase H2O. In 
particular, we find PBE GqWq ( 6.0 eV) ~ PBE DFT < HSE 06 GqWq 
^ PBE scQPGIVl < HSE 06 DFT ( 7.3 eV). Figureshows that, 
as was the case for H2O in gas phase, the IP of the ll 20 @Ticus 
interface across the various methods is ordered according to the 
method’s description of the screening, . 

As discussed above, for hybrid xc-functionals such as HSE, the 
effective screening is determined by the fraction of exact exchange 
a included. Essentially, a plays the role of the effective screening 
within the method, s '^. Although HSE 06 incorporates less screen¬ 
ing (^oo ^ 4 ) than experiment for rutile Ti02 (£™^ ^ 7 . 6 ),^ the 
HSE 06 IP for the inter face i s in agreement with the experimental 
estimate of /P 7.1 

If one performs GqWq based on HSE 06 , a stronger screening is 
applied, i.e., £00 5 . 7 , yielding a lower IP for the interface. In fact, 

as indicated by the red arrow in Figure a similar IP to HSE 06 
Go Wo should be obtained from HSE DFT by setting the fraction 
of exact exchange to the inverse dielectric constant of bulk Ti02, 
i.e., a = 1 /e™^. Adjusting a to the measured inverse dielectric 
constant has been previously found to give improved band gaps.® 
From PBE scQPGWl, one obtains an IP consistent with that of 
HSE 06 Go Wo . This is because we find the screening in scQPGW 
decreases from PBE RPA with each self-consistent cycle. Essen¬ 
tially, the final screening incorporated in scQPGWl is similar to 
thatofHSE 06 RPA. 

As shown in Figure [TT] PBE Go Wo gives an IP slightly lower 
than PBE DFT for the interface, while the PBE Go Wo CBM is 
shifted up by about 2 eV. This is surprising, since PBE DFT already 
yields a CBM level alignment for the interface in excellent agree¬ 
ment with experiment. This is partially due to PBE RPA’s over¬ 
estimation of the screening of Ti02 (^oo ~ 8 . 3 ). Although HSE 06 
Go Wo has a weaker screening than PBE GqWq, the resulting abso¬ 
lute alignment of the CBM is quite similar. If instead, the self en¬ 
ergy corrections are applied self-consistently via PBE scQPGWl, 
the absolute alignment of the CBM is significantly lower, but still 
greater than that of PBE DFT or HSE 06 DFT. This is again related 
to decreases in the dielectric constant with each self-consistent cy¬ 
cle. For this reason, scQPGWl tends to provide reasonable band 
gaps for TiO2(110) interfaces. Overall, we observe an ordering in 
increasing band gap of PBE DFT < HSE 06 DFT < PBE scQPGWl 
< PBE Go Wo ^ HSE 06 GqWq, with HSE 06 DFT providing the best 
absolute alignment of the CBM and VBM for the H20@TicM^ in¬ 
terface. 

In Figure we show that a similar correlation between IP 
and the method’s description of screening holds for clean and 
hydroxylated h-TiO'j(llO). Specifically, we consider clean stoi¬ 
chiometric TiO2(110),^and dissociated H20@0^^^ on defective 
Ti02-i/4(l 10 ) with ViML of Overall, IP ^ /Ppbe + 5 . 9 a for all 


systems considered. We again find that the IP of PBE Go Wo ~ PBE 
DFT, HSE 06 GqWq ~ HSE(Qf“' = pbe scGWl, and PBE 

scQPGW HSE 06 scQPGW ~ HSE 06 DPT. 

HSE 06 DFT provides the most accurate description of the IP of 
the clean and ll20@Ticus covered stoichiometric TiO2(110) sur¬ 
faces. Although the HSE 06 DFT IP for H20@0^^‘^ is significantly 
lower than the one measured for /z-TiO2(110), in both cases, the 
IP is shifted to lower energies relative to the clean stoichiometric 
surface. Differences in the magnitude of the shifts are probably 
due to the differences in defect coverage between the experiment 
( 6 - 9 %)® and the calculation ( 50 %). 

The similarty between HSE 06 DFT and scQPGW based on either 
PBE or HSE 06 for the clean TiO2(110) surface,® points to a sim¬ 
ilar screening from these two techniques. This also demonstrates 
the starting point independence of the scQPGW technique. 

To summarize, although scQPGW provides accurat e IPs, the 
band gap is greatly overestimated, as reported previously.* —*—*—* 
While scQPGWl provides a more accurate band gap, it achieves 
only a qualitative description of the IP. HSE 06 achieves a quan¬ 
titative description of both the IP and band gap, but provides a 
poor descripti on of the molecular level alignment relative to the 
YBM ISmEll the highest occupied H2O levels are 

significantly hybridized with the substrate, this is not a major draw¬ 
back in this case. In general, for TiO2(110), a more effective strat¬ 
egy is to combine the calculated IP from HSE 06 with the occupied 
interfacial levels’ alignment from GqWo or scQPGWl. 

4. CONCLUSIONS 

The level alignment prior to photo-irradiation is an important piece 
of the puzzle needed to get a complete atomistic picture of pho- 
tocatalytic processes. Here we have shown that the complex UPS 
spectra for the H20-Ti02 interface may be disentangled using QP 
Go Wo PDOS. We have firmly established the robustness of the QP 
Go Wo H2O PDOS by: ( 1 ) demonstrating its xc-functional (PBE, 
EDA, vdW-DF, and HSE 06 ) independence, ( 2 ) comparing to self- 
consistent QP GW techniques (scQPGWl), and ( 3 ) considering its 
dependence on surface coverage and dissociation. Altogether, these 
calculations provide an accurate interpretation of the complex UPS 
and MIES experiments®®^ for the H2O-TiO2(110) interface, 
and provide accurate estimates of the highest H2O occupied lev¬ 
els’ alignment relative to the VBM. 

Our results provide two important pieces of the puzzle: ( 1 ) 
the molecular structure of the photocatalytic interface and (2) the 
molecular alignment of the doubly occupied levels near the VBM 
responsible for hole trapping prior to irradiation. To complete the 
picture, the molecular structure and level alignment in the pres¬ 
ence of the photo-generated hole is also needed. Previous DFT 
studies using the hybrid HSE xc-functional have found a hole can 
be trapped at surface O 2 p;r levels of and HO-Ticw,? sites. ^ 
However, the screening of such localized levels may not be well 

described by HSE, which tends to underbind localized interfacial 

I 23 I 

levels.* — 'This underbinding is corrected upon inclusion of many- 
body effects via QP GoWq.® Having demonstrated the capability 
of GoWo for the description of level alignment prior to irradiation, 
this work points the way forward via future QP GoWq studies of 
level alignment for trapped hole levels. 
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